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Abstract 

We investigate the electromagnetic properties of the deuteron such as the charge and 
magnetic form factors by solving the Bethe-Salpeter equation (BSE) with the separable 
ansatz. In solving the deuteron bound state solution to the BSE, we include negative 
energy components of P-wave in addition to the 3 Si and 3 Di states of positive energy, 
employing a rank IV separable ansatz. We found that the inclusion of the negative 
energy components improves systematically the electromagnetic properties which are not 
described in the conventional non-relativistic impulse approximation. 

1 Introduction 

Traditionally, nucleon-nucleon scattering and deuteron properties such as electromagnetic 
properties are investigated in the non-relativistic framework. For example, Gari and 
Hyuga investigated electromagnetic properties of the deuteron up to momentum transfer 
2.4 GeV 2 [TJ [2]. In order to explain the experimental data, they needed to introduce 
meson-exchange currents in addition to the impulse processes. In 90's, Tamura et al., 
Blunden et al. and Wringa et al. j3jlU[5] a l so investigated the electromagnetic properties of 
the deuteron in a non-relativistic framework. In the non-relativistic framework, electron- 
deuteron scattering is expressed as shown in Fig. [H where (a) is the impulse process 
and (b) is the contributions from the exchange currents. However, all non-relativistic 
treatments which are based on the impulse approximation with relativistic corrections 
and meson exchange currents are sensitive to the model of the meson exchange currents. 

In principle nucleons are described by the Dirac equation, which is the relativistic 
equation. In the relativistic framework, we expect the following advantages: 

• Relativistic kinematics are automatically included. 

• Some exchange currents are included through the Z-graph process. 

• Dynamics originated from the relativistic effect such as the LS-force, spin-spin in- 
teraction and the negative energy components are naturally included. 

In the relativistic framework, electron-deuteron scattering can be expressed as shown 
in Fig. [2j In the relativistic impulse approximation (RIA), each nucleon line has both 
negative and positive energy states, in contrast with the non-relativistic approximation 
where each nucleon line contains only positive energy states. The right hand side of Fig. [2] 
shows the decomposition of the RIA diagram into several diagrams in the non-relativistic 
reduction, where the diagrams are depicted in the chronological order. One of those 
diagrams expresses the Z-graph whose nucleon line goes back as corresponding to the 
negative energy states. Therefore, we can interpret that the Z-graph part corresponds to 
the exchange current in the non-relativistic approach. It implies that a part of dynamics 
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of the inclusion of the negative energy components can be interpreted by the exchange 
currents. 

In order to see the effects of the relativistic framework, we treat the Bethe-Salpeter 
equation in the fully relativistic manner [6j. In 1989, Rupp and Tjon solved the Bethe- 
Salpeter equation (BSE) by the covariant Graz II potential. It is a simple covariantization 
of the Graz II potential which was built by Mathelisch, Plessas and Schweiger [3 [8]. Rupp 
and Tjon succeeded in reproducing 1 Sq and 3 £i- 3 .Di phase shifts up to 600 MeV in their 
relativistic framework. However, they did not include negative energy components. For 
the investigation based on a fully relativistic framework, we should include the negative 
energy states appropriately. 

For the relativistic description of the deuteron, we prepare a two-nucleon set including 
negative energy states of P-wave in addition to the 3 S± and 3 D\ states of positive energy. 
We solve the BSE by using the separable ansatz including the negative energy states, and 
determine the parameters to reproduce deuteron form factors. 

The contents of this paper are as follows. In section 2, we give a general formalism of 
the Bethe-Salpeter equation for the deuteron, and provide a solution using the separable 
ansatz. We discuss how to determine the parameters of separable potential in some 
detail. In section 3, we show the relativistic kinematics for the elastic electron-deuteron 
scattering. In section 4, results and discussions are presented. The final section is devoted 
to the conclusions of the present work. 




(a)NRIA (b)Exchange current 



Figure 1: Feynman diagrams for non-relativistic impulse approximation (NRIA). 




NRIA Exchange current 



Figure 2: Feynman diagrams for relativistic impulse approximation (RIA) and its 
non-relativistic decomposition. 



2 BS approach with a separable interaction 
2.1 The BS equation 

Let us start with the BSE for the NN T-matrix: 

+i J 4^V aP , eX (P,p', fc)5g)(P/2 + k)£$(P/2 - k)T mSj (P, k,p) , (2.1) 
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where the Greek letters express spinor indices. The T-matrix and the interaction kernel 
V are the functions of the total momentum P, and the relative momentum p and p', of 
the initial and final states. is the one nucleon free propagator. The bound state 

corresponds to a pole in the T-matrix at P 2 = M B , where Mb is the mass of the bound 
state: 



Ta.ppfsi.Pi k ,k) — ^ p 2 ^ ^ 2 + Ra/3,yS (-P> k j k) 



(2.2) 



Here Y a p is the vertex function of BSE, and R a 0^& is regular at P 2 = M B . We can 
express the BS amplitude by the vertex function as : 



P 



P 



^(P, k) = S$(^ + k ) S $C-2 ~ k)T ^ k) ' 
and we obtain the equation for the BS amplitude from Eqs. (|2.1j) . (|2.2p and ()2.3f) : 

,P „„m,P 



(2.3) 



$ a(3 (P,k)=iS£j(- + k)S$(--k) J —V meX (P,k,k")*ex(P,k"). (2.4) 

2.2 Solutions to the BSE 



After partial wave decompositions (see [E]), the BS equation for the T-matrix in the 
center-of-mass frame is given by 



T a p(p' , \p'\,Po, \p\',s) = V a p(p' , \p'\,p Q , \p\;s) + 



(2.5) 



% 

2^2 



j dk k 2 d\k\ Varf(p'o, \p'\,k , \k\;s) S^s(k , \k\;s) T S)3 (k , \k\,p , \p\;s) . 

■yS 



Here the Greek indices represent quantum numbers (JLSp), and the summation takes 
over all partial waves. When we include both positive and negative energy states, we have 
eight states, namely, 3 S'J t ", 3 T> 1 f , 3 ST, 3 T>7, 1 Pf, 1 P°, 3 Pf and 3 P°, which are labeled as 
1, . . . , 8. Details of the structure of these amplitudes are discussed in [16]. Therefore T, 
V and S are 8x8 matrices where n-th row and column corresponds to the n-th state. In 
this eight dimensional basis, the propagator S a p is expressed by 
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where two nucleon propagators S a are given as 

1 



S 4 



(s/s/2 -e k + k + ie) (y/s/2 - e k - k + ie) 
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S- = 



1 



1 



{y/s/2 + e k + k - ie) {y/s/2 + e k - k - ie) ' 



Q — 

{y/s/2 - e k + ko + ie) {y/s/2 + e k - k - ie) ' 

1 1 

{y/s/2 + e k + k - ie) {y/s/2 - e k - k + ie) ' 

_ _ _ S-\ — + S — |- 

''(: "■'(!(! ^OO 2 ' 

g _ g 

So = S eo = S oe = - . (2-7) 
Now let us introduce the separable ansatz of rank N in the following manner: 

N 

V a f3(Po, \p'\,Po, \p\;s) = E (Po, \p'\)9j (Po,\p\)> \j = Xji, (2.8) 

i,i=l 

where we assume that Ajj is symmetric under the interchange of Then the T-matrix 
is also given in a separable form as 

N 

T af} {p' ,\p'\,p ,\p\;s) = E n ] gt\p^\p'\)9f\pM), (2-9) 

and 

{T- 1 {s)) ij = {\~ 1 ) ij -H ij {s). (2.10) 

Here Hij{s) is defined by 

^i( s ) = ^ E [ dk ° k2 d \ k \ V( fc o, |fe|; s) gf LSp) {ko, \k\)g( JLSp '\k , \k\). (2.11) 

/7T LSpp' J 

Then, the radial part of the BS amplitude can be written as 

4>jls p {po, \p\) = J2s pp >{p ,\p\) g.jLs P {po,\p\) 
p' 

N 

= EE S ppl {pM;s)X l] g\ JLSp '\p 0l \ P \y 3 {s)- ( 2 - 12 ) 
p' ».j=i 

where q(s) satisfy the following equation: 

2V 

C-(8) - E H ik {s)X kjCj {s) = 0, (2.13) 



Using (j)jLSp{Po, \p\), w e can obtain the general form of the BS amplitude (see P33]). 

Because 3 Pi spin- angular part has a complex form (see [E]), we include *Pf and 
1 Pf but ignore 3 P° and 3 Pf among negative energy components for simplicity in actual 
calculation. We also ignore 3 5 , jf, 3 -D 1 " states, because we expect that the contribution from 
those two states is second order in the non-relativistic expansion (in powers of velocity) 
and is small. 
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2.3 Rank IV separable ansatz 

For a relativistic description of the deuteron when we include both the positive and nega- 
tive energy components in the BS amplitude, we extend the covariant Graz-II interaction 
with including P-wave components. For this purpose, we adopt a rank IV separable 
ansatz. The Graz-II interaction is a rank III separable potential suggested by Mathelitsch 
et al. in the non-relativistic framework [7j. Rupp made a covariant Graz-II potential 
by simple covariantization [8]. In the covariant Graz-II interaction, they considered only 
positive energy states, 3 Sf and 3 D^. Here in this work we build a rank IV separable po- 
tential (N = 4) by improving the rank III covariant Graz-II interaction with the inclusion 
of the negative energy states. We use the following forms for the functions gf, 

9i (po,\P\) ~ 



3 s+ , h (Po-P 2 ) 



[pi -p 2 - 0f 2 + ief 



93 (P °' IPl)= (pl-pi-^+ieM-pi-^ + ieT (2 - U) 



lp if i i\ p 
9a (po,\p\) 



9a 1 (Po, \p\) = 73 



(p§ -P 2 -f3l + it) 2 ' 
Po \p\ 



' 771 (Po ~ P 2 ~ Pi + 7e ) 2 ' 

3g+ 3p) + lpe,° 

93,4 (P0» |p|) = 01,2,4 (PO, \P\) = 9l,2,3 (P0> \p\) = 0, 
01,2*3,4 (PO, \P\) = 01,2,3,4 (PO, \p\) = 01,2*3,4 (PO, \p\) = . 

For 3 Sf and 3 Df, we use the same type of function as the covariant Graz-II Interaction. 
What is new here is the inclusion of functions for the 1 P^'°-waves. The last line reflects 
our assumption that we ignore 3 P° , 3 Pf, 3 5'~ and 3 -tq . 

Having the ansatz Eq. (|2.14p . the solution Eq. (|2.12p to the BSE can be written as 

4>s s +{Po,\p\) = (ciAn + c 2 Ai2 + C3A13 + c A Xu)S + g 1 1 (p , \p\) + (2.15) 
(C1A12 + c 2 A 2 2 + C3A23 + c 4 A 2 4)5'+92 1 (Po, \p\), 

3 D + 

4>3 D +(po, \p\) = (ciAi 3 + C2A23 + C3A33 + c 4 A 34 )5' + c/3 1 (po,\p\), 

1 pe 1 po 

fapf(Po,\p\) = (ciAu + c 2 A 2 4 + C3A34 + c 4 A 44 )(5' e 5 4 1 (po, \p\) + Sog^ 1 (p , \p\)), 

1 po 1 pe 

4>ip°(Po,\p\) = (ciAu + c 2 A 24 + c 3 A 34 + c 4 X 4A )(S e g A 1 (p , \p\) + S g 4 1 (p , \p\)) , 

03 S -(PQ, \P\) = <P3 D -(P0, \P\) = 4>3po(p , \p\) = (f)3pe(po, |f>|) = , 

where we have used the facts that Xij = Xji, and that (pipe is even and that (ftipo is odd 
under po — > —po. 

2.4 Computations with negative energy components 

When we use the rank IV separable potential, Eq. (|2.13p takes the following form 

3 

ci — ^ HikXkjCj = ca(H n Xia + H 12X24) , 
k,j=l 
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3 

C2 — ^2 Hlk^kjCj = C4(H2\Xl4 + H22X24) , 
k,j=l 
3 

C3 — X! Hsk^kjCj = 04(1133X34) , 

k,j=l 
3 

C4- ^2 H ^kXkjCj = C 4 (if4 4 A44) . ( 2 -16) 
k,j=l 

In these equations, Xj4, H^jij = 1,2,3,4) and C4 appear due to the inclusion of P-f and 
1 P° -states. Since the fourth component in the rank IV ansatz is introduced for the first 
time in this work, we would like to discuss in some detail the quantities carrying the 
fourth index. First, we note that i?4j(j = 1, 2, 3) = due to the choice of the g- functions 
as given in Eqs. (|2.14p . On the other hand, #44 can take a finite value. We would like to 
discuss H44 in detail. 

The evaluation of #44 contains a /co-integral, an d the result is affected by the location 
of poles of the ^-functions and propagator S a as defined by Eq. (|2.7p . For deuteron 
problems, we take the location of poles of g 1 and g p ° in the same side of the complex 
fco-plane as g Sl . As we will discuss in detail later, P-wave part influence on the form 
factors and tensor polarizations through #44, Cj(j = 1, 2, 3, 4) and Xjj(j = 1, 2, 3, 4). This 
is an important contribution in order to improve the agreement with the experimental 
data through the negative energy components in bound state problems. 

However, if we adopt the same prescription for scattering problems, a finite #44 causes 
a serious problem, which change the phase shift drastically, even if the rate of -P-wave is 
very small. Therefore, we try to solve this problem in the following way. For H44, we 
have the following expression: 

H44 = i / dk ° fc2 d \ k \ ( S e + {g\ P? ?) + 2S [g\ P * g/ ? ^j 

= ^2/ dk k 2 d\k\ (s + -(g + -) 2 + S- + (g- + ) 2 ) , (2.17) 

where 

Ipe 1 po 

9 =g 4 1 + 94 1 > 

1 pe 1 po 
h 1 1 

9 =94 ~94 ■ 

At first sight, Eq. (|2.17p takes a finite value. However, by locating the poles of g^ in 
the upper side of the complex /co-plane as , and the poles of g h lower as S h (see 
Fig- El), Eq. (|2,17p becomes zero. This prescription to make the vanishing contribution 
for H44 turns out to be crucially important in order to reproduce the experimental phase 
shifts (see Fig. Physically, this means that we neglect the direct contribution from 
negative energy components for scattering problems. For the present separable potential 
of g-functions, we have realized this by locating the pole in the appropriate manner. As we 
discuss in Appendix [Bl the similar situation can be realized for static separable potential, 
which justifies the present prescription. 

By assuming different pole location of the ^-function for the bound state and scat- 
tering problems, we are able to reproduce reasonable results for both quantities. The 
physical ground of this condition is related to the fact that a finite value of H44 allows 
appearance of negative energy states in intermediate states. For a bound state, negative 
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Figure 3: Location of poles of S+_ and S-+ 



energy components are naturally contained in the relativistic formalism, and therefore, 
the mixture of negative energy state is allowed. In contrast, for scattering problems such 
mixture is not allowed, since scattering states are expressed as definite energy states as in 
the asymptotic states. Therefore, H44 must vanish for scattering problems. In this way, 
we will be able to obtain good results for both scattering and bound state (deuteron) 
problems simultaneously. 

In actual calculations, we express the Aj4 parameters in terms of a single 114 parameters 

by 

A14 = — \/\iiU4, A24 = VA22U4, A34 = V/A33U4, A44 = u\ . (2-18) 

We try to find good parameter sets to reproduce the form factors, tensor polarizations 
and phase shifts. 




Figure 4: Phase shifts of 3 Si (left) and 3 P>i (right) channels. Data are taken 
from SAID program (http:/ /gwdac.phys.gwu.edu). Calculational results (solid 
line) are performed with the D-wave ratio 4% and with the negative energy 
P-wave ratio 0,-1 and —3% (for negative values of the P-wave ratio, see the 
discussion in Appendix C). Results with three P-wave ratios are almost the 
same within the resolution of the figures. 
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3 Elastic Electron-Deuteron Scattering 
3.1 Relativistic Kinematics 

The differential cross section for unpolarized elastic electron-deuteron scattering in the 
one-photon-exchange approximation is given by 

% = (^-\ U(s 2 )+B( 9 2 )tan 2 !l. (3.1) 
Here ( 4S- ) is the Mott cross section given by 

V dn o J Mott to J 

da \ a 2 cos 2 9 e /2 



dtt' e J Mott AE 2 (1 + 2£ e /M sin 4 e /2) ' 



(3.2) 



where 6 e is the electron scattering angle, M the deuteron mass and E e the incident electron 
energy. The functions A(q 2 ) and B(q 2 ) are the deuteron structure functions which can be 
related to the form factors by 

A( q 2 ) = F 2 ( q 2 ) + S - V 2 F 2 l ( q 2 ) + 2 - V F 2 1 ( q 2 ), 

B(q 2 ) = ^r ] (l+r ] )F 2 l (q 2 ), (3.3) 

where r\ = —q 2 /AM 2 = Q 2 /AM 2 . The electric Fc(q 2 ), the magnetic Fu(q 2 ) and the 
quadrupole FQ(q 2 ) form factors are normalized as 

F C (0) = 1, F m (0)=^d — , F Q (0) = M 2 Q D , (3.4) 

m 

where m is the mass of the nucleon, jUd the magnetic moment and Qd the quadrupole 
moment of the deuteron. The tensor polarization components of the final deuteron can 
be written through the deuteron form factors as follows: 

t 20 [A + man 2 |] = -l=[^r,F c F Q + ^ V 2 F 2 + ^(1 + 2(1 + rj) tan 2 ^)F 2 ], 
«21 [A + Btzn 2 ^] = -j=v(ri + r/ 2 sin 2 9 -^) l / 2 F u F Q sec |, (3.5) 



M + ^la.i 2 y] =-^=??F M . 



We can obtain Eq. (|3.ip from the following amplitude 

Af fi = ie 2 u m ,(l')^u m (l) \ (D'M'\J^\DM), (3.6) 

where u m (l) is the free electron spinor with 4-momentum I and the spin projection m, 
q = I — V = P' — P the 4-momentum transfer and P (P 1 ) the initial (final) deuteron 
momentum. \DM) is the deuteron state with an angular momentum projection M and 
is the electromagnetic current operator. 

The deuteron current matrix element is parameterized in the following way 



{D'M'\J»\DM)= eCaM'iP') 



(P' + P^F^-^Fziq 2 ) 



&M(P) , (3-7) 



S 




Figure 5: A schematic diagram for the elec- 
tromagnetic process of the deuteron. 

where £m (P) and ^* M , (P') are the polarization 4- vectors of the initial and final deuteron. 
Form factors Fi^q 2 ), G\{q 2 ) are related to Fc(q 2 ), Fq(q 2 ) and FM_(q 2 ) by the equations 

^0 = ^1 + \v\Fi + (1 + V)F 2 ~ Gi] ,F Q = F 1 + (l + rj)F 2 -G u F M = G l . (3.8) 
The normalization condition for the deuteron current matrix element is given by 

lim (D'M'\J^\DM) = 2eP IM S MM > . 

Equivalently, in terms of the BS amplitude: 

/ (0^ f(P ' k) JP, { S(1) ( P / 2 + k )S (1) ( P / 2 ~ *)} W *) = -2P, . (3.9) 

In order to calculate the deuteron form factors, let us choose the laboratory frame 
(deuteron at rest). In the laboratory frame the relevant momentum variables take the 
following form (the z-axis is along the photon momentum): 



P = (M, 0), P' = P + q = (M (1 + 2??), 0, 0, 2MVWTT^), 

q = (2Mr j ,0,0,2M y /Tjy/T+v), (3-10) 

£ M=+1 (P) = Cm=+i (p') = - -j= (0, 1, i , 0) , 
£ M =-i(P) = Zm=-i(P') = ^=(0,1, -i,0), 



tiM=o(P) = (0,0,0,1), £ M =o(P') = (2^71+^,0,0,1 + 27?). (3.11) 
From Eqs. (HTTPD . (|3TT]) and Eq. $Tfl), we obtain: 

(M'\Jo\M) = 2Me (1 + t?) {f 1 5 mm > + 2r?[Fi + (1 + r?)F 2 - Gi]Wmo}, 
(M'\J X \M) = ^£ VV V^+V G x {S M 'M+i ~ S M 'M-i}- (3.12) 

3.2 Deuteron currents and form factors 

In the relativistic impulse approximation, a deuteron current matrix element can be writ- 
ten as 

(D'MV^DM) = iej ^Tr[^M'(P , ,p')T^ +n Hq)^M(P,p)iS^ T (q 2 ))- 1 

r,(q) = l,F x {q 2 )-l^p±F 2 {q 2 l (3.13) 
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where (P, p) is the BS amplitude of the deuteron, P' = P + q and p' = p + q/2, where 
q is the momentum transfer. T^(q) is the vertex of 'yNN interaction, F\ and F2 are the 
Dirac and Pauli form factors of the nucleon, and the sum over the proton and neutron is 
taken. A schematic diagrams is shown in Fig. [5j 

The resulting expressions for the deuteron current matrix element can be written as 



(D'M'\J^ IA \DM)=l^ M (q 2 ) Fife 2 ) + 2jft*V) F 2 (q 



rM'M/JZ 



M'M 1 



(3.14) 



rM'M(„2 



1,2 M 
L'S'p',LSp 



J dko \k\ 2 d\k\ d(cosl 



fco,|fc|) I 



L',L 

1,2 M'A^ ^ 



(fco, |fc|, cos g 2 



1,2 M'M At 



(fco, |fc|, cos 6, q 2 ) is obtained by taking the trace in the 7 



where the function / 

matrix space of Eq. (|3. 13j) and the substitution of the scalar products into Eq. f)3. 13j) . 

In Eq. (|3.14j) . the radial part of the BS amplitude for the final state deuteron (f>ws'g'(k' Q , \k\ 
depends on the momentum variable k' in the laboratory frame. The momenta of the ini- 
tial deuteron (P) and of the final deuteron (P') in the laboratory frame are related by 
the Lorentz transformation: 



p> = £P = £(M,0), k' = £k, 
where the Lorentz transformation matrix £ is of the form: 

( I + 277 2y/rjy/T+Tj ^ 
10 

1 

1 + 2t? / 



(3.15) 



(3.16) 



V 2^7] v / m? 

The components of the 4- vector of the final state k' are expressed by k' = (k' , k' x ,k' y , k' z 
and |Jfe'| = Jk'J 2 + k' 2 + k' z 2 . Using Eq. (pTTD]) . (13T51) and (ETI51) we obtain 



A, 



(1 + 2 V )k - 20^1 + i]k z - M V , 
kxi ky = ky, 



(1 + 2rj)k z - 2^1 + 7/fco + M07vT + r7, 



(3.17) 



where ko,k x ,k y , k z are the components of the 4- vector k of the initial state. 

In order to calculate the deuteron form factors, we need to know three matrix elements 
with different total angular momentum projections and current components, for example, 
(0|Jo|0), (1|Jq|1) and (1| J x \0). Finally we obtain the following equations. 



Fc 
F M 
Fn 



1 (p'M 1 = 0| J Q \PM = 0) + 2(P'M' = +1| J \PM = +1) 



2M 
1 (P'M' 



3(1 + 7?) 

+1\J X \PM = 0) 



My/2 ^qy/TTv 
1 (P'M' = 0| J \PM = 0) - (P'M' = +1| J \PM = +1) 



(3.18) 



2M r/vTTr? 
We can obtain tensor polarization components by inserting these equations into Eqs. (|3.6|) 
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4 Results 



We have calculated the deuteron form factors and tensor polarizations. Our purpose 
is to study effects of negative energy components. Therefore, first we show results of 
the impulse approximation without negative energy components. As shown in Figs. [6l 
the results of this impulse approximation does not reproduce experimental data. Then 
we attempt to improve the agreement by changing nucleon form factors and D-wave 
probability. However, as shown in Figs. El the two effects do not improve the agreement. 

4.1 Dependence on the intrinsic nucleon form factor 

We consider the typical three types of nucleon form factors: the dipole fit, the vector meson 
dominance model (VMDM) and the relativistic harmonic oscillator model (RHOM). The 
details of these form factors are presented in [16]. The three form factors reproduce the 
proton charge radius equally well. However, they differ in the momentum dependence 
and in neutron form factors. Therefore, our interest is how these differences affect the 
deuteron form factors. In order to see the effect of the different form factors, it is sufficient 
to show the results which are calculated including only positive energy states, 3 Sf and 
3 -D^, and with the ratio of the D-wave being fixed at 4%. (The evaluation of probabilities 
of the BS amplitudes is discussed in Appendix [A]) For this probability, we use the same 
parameter set as used by Rupp and Tjon [8] for Pu = 4%. In our case, we have additional 
parameters Aj4 (i = 1, 2, 3, 4) and 73, which are set equal zero when the P-wave is ignored. 

Results are shown in Fig.[6}-(a) for the charge form factors of the deuteron, where solid 
line is calculated using the dipole form factor, the dashed line with the VMDM, and the 
dotted line with the RHOM, which are compared with experimental data [TT]. From this 
comparison, we verify that the deuteron form factors are not very sensitive to the nucleon 
form factors. This result is understandable since the major structure of the deuteron form 
factor is determined by the loosely bound nature of the wave function which makes the 
form factor fall off very rapidly as the momentum transfer Q 2 is increased. As compared 
to the wide spreading structure of the deuteron wave function, nucleons are regarded as 
small objects. Although we have shown the result of the charge form factor only, other 
quantities such as magnetic and quadrupole form factors are also not sensitive to the 
models of the nucleon form factor. 

4.2 The role of the D-wscve 

We investigate the role of the D-wave in the charge form factor. The D-wave ratio is 
usually estimated as 4 ~ 6%. Therefore, we have calculated form factors with the D-wave 
ratio at 4,5 and 6%. We use the parameter sets of Rupp and Tjon for Pp = 4, 5 and 
6%. As in the previous subsection, parameters Xa (i = 1, 2, 3, 4) and 73 are set equal zero. 
Results are shown in Fig. [H-(b). The solid line is the result obtained with the D-wave 
ratio 4%, the dashed line 5%, the dotted line 6%. From this, it is not easy to make a 
significant improvement as to reproduce the experimental data with a reasonable range 
of D-wave ratio. 

4.3 The role of negative energy P- waves 

In the previous subsections, it was shown that nucleon form factors and the D-wave ratio 
did not improve the agreement with the experimental data. Therefore, in this subsection, 
we investigate the effect of the negative energy P-waves on the form factors and tensor 
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polarizations. We performed calculations by including the two P-wave components of Pf 
and 1 P° as explained in the section 2.2. We adopt the nucleon form factor of dipole fit 
and the D-wave ratio 4% 

The results are shown in Figs. [JJ Fig. [7]-(a), (b), (c) show charge, magnetic and 
quadrupole form factors respectively, and Fig. [TJ-(d), (e), (f) the i20> *2l and £22 tensor 
polarizations, respectively. 

The dotted lines represent the result when the the P-wave ratio is 0%, the dashed 
lines —1% and the solid lines —3%, except for Fig. H-(b). For Fig. [JJ-(b), the dotted, 
dashed and solid lines represent the results with P-wave ratio is 0%, —0.5% and — 1%, 
respectively. The negative values are called pseudo probabilities which are related to the 
baryon charge [2]. Details are discussed in Appendix C. The ratio of negative energy 
states is controlled by the parameter which is chosen as M4 = —3.5, —5.0 and —8.5 
for the P-wave ratio —0.5%, —1% and —3%, respectively. The corresponding g function 
parameters are chosen as = 0.481 GeV and 73 = —15.0. 73 and are determined 
as follows. We try to reproduce the location of the dip (sign change) of Fq, firstly by 
setting 73 = —15.0. (3% is then determined so that H^s = M 2 ) = 0. The rest of the 
parameters are taken from Rupp and Tjon Ref. [8] of D-w&ve ratio 4%. The data of the 
charge and quadrupole form factors are taken from Ref. [11], those of the magnetic form 
factor Ref. |13j . and those of the tensor polarizations Refs. |11|, I12j. 

We find that better agreement with experimental data is achieved when we include 
a finite amount of P-wave amplitudes of the negative energy components. For charge 
and magnetic form factors, the location of the dip and the Q 2 dependence especially at 
high momentum region Q 2 > 1 GeV are significantly improved by including the negative 
energy component. The amount of the negative energy ratio differs, however. The best 
agreement is achieved with —3% for the charge form factor, while —0.5% for the magnetic 
form factor. No significant change was found as for the quadrupole form factor. As for 
the tensor polarization, the agreement for £20 is remarkably improved when the negative 
energy ratio is —3%. The agreement of the other two components is also reasonably 
improved with a finite amount of negative energy components. The necessity of different 
amount of negative energy components in order to reproduce different form factors might 
be a consequence of our incomplete treatment of the negative energy components, where 
we have ignored negative energy components of 3 Pf , 3 Pf , s Si and . The inclusion 
of the full set of components will be a future work. 

5 Conclusion 

In this paper, we have investigated the effect of negative energy components on two- 
nucleon systems. For that purposes, we have solved the Bethe-Salpeter equation, including 
the 1 Pf, 1 Pf states, both for scattering states and the bound state. We have found 
the parameter sets which reproduce electromagnetic properties and the phase shifts of 
the deuteron. On one hand, the inclusion of the negative energy components improves 
systematically bound state properties such as the form factors and tensor polarizations 
in the relativistic impulse approximation. On the other hand, using the same parameter 
set we have seen that the negative energy components affect very little on scattering 
properties. 

Although we have not discussed in detail, we have compared our results with the 
previous study based on the non-relativistic formalism. In the non-relativistic method, 
the results depend very much on the approximation, the way how exchange currents 
are included, and so on. In the present relativistic method, although our scheme is a 
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Figure 6: The Charge form factor of the deuteron as function of Q 2 . The 
results of the left panel (a) are calculated with the D-wave ratio 4% and with 
three different form factors, Dipole, VMDM and RHOM (see text). The results 
of the right panel (b) are calculated with the nucleon form factors Dipole and 
with the D-wave ratio 4%, 5% and 6%. For both cases, the negative energy 
P-wave ratio is fixed to be zero. 



simple impulse approximation, the inclusion of negative energy components systematically 
improved the agreement with data. In this respect, it would be of great interest to 
investigate the non-relativistic correspondence of the present relativistic method. 

Although the present approach provides formally a systematic method to study two 
nucleon systems in a relativistic way, actual computations are rather complicated due to 
the increase in the number of components of the BS amplitudes and the analytic structure 
in the complex energy plane. For the reasons associated with these practical problems, 
we had to ignore mathematical rigor where only part of negative energy components were 
taken care of. Nevertheless, our attempt to include for the first time negative energy 
components have shown that it indeed improves the agreement between the theory and 
experimental data in a systematic manner. It will be of great interest to further investigate 
the role of negative energy components, which is equivalent to the antiparticle degrees of 
freedom. Such dynamics is expected to become important for more strongly interacting 
systems such as quark systems for hadrons. 
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Appendix A Probabilities of the BS amplitude 

In this appendix we discuss how probability is defined for the BS amplitude. From the 
normalization condition Eq. (|3.9p when /i = 0, we can obtain the normalization condition 
for the radial part as follows, 



dkr, 



d\k\\k\ 2 r\k ,\k\ 







S(k , \k\;s) 



r(*o,|fc| 



-2M, (A.l) 



where the integration over angles has been done and r(/co, |fe|) is the radial part of the 
vertex function. For the deuteron of J = 1, S is the two nucleon propagator define by 
Eqs. (|2.6p and (|2.7p which is given as a matrix form for eight channels. 
Substituting the definition of the radial part of the BS amplitude 



rfo, \k\) = s-^ih, \k\) 

obtained from Eq. (|2.3|) into Eq. (|A.ip . we find 
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d|fc||fc|V(fco,|fc|)[S ,_1 ; 
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J dk o J d\k\\k\ 2 <j)* a (ko, \k\)u a p(j)p(k , \k\) 



= -2M. 

In the second line we have used eight components explicitly, 
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f 35 +(*o,|fc|)^ 
4> SD +{k , \k\) 
(t> 3s -(k ,\k\) 
(/> 3D -(k ,\k\) 

01pe(ko, |fc|) 

<j)apo(ko, \k\) 
<fnpo(k , |fc|) 
V <t>3pe(k , \k\) J 



4>\k , \k\) = {(p* 3s +{k , \k\),<fe D +(ko, 1*1), (f>t s - (k , |*|), 4>t D - (*d, 1*1), 

0*pe(/SO, |fc|), ^3 po(fe , |fe|),^*po(*0, |*|),^3pe(A;0, |*|)) • 

From these equations, we can define the probability of finding the a-state in the total 
deuteron state by the following expression, 

P a = ^J dk J d\k\\k\ 2 u; a \cf> a (k ,\k\)\ 2 , (A.4) 

where N is determined by 

£ P a = 1 (a = 3 S+ 3 £>+ 3 Sf, 3 Df, 1 Pf, 1 Pf, 3 Pf, 3 Pf) . (A.5) 

a 

P a is often called as "pseudo probability" which is normalized by the baryon charge [14j . 
When a = 3 Sf, 3 Pf, 1 Pf , Pf, 3 Pf , Pf , the sign of P a is minus in accordance with the 
sign of uiaaip- = 1, - - - , 8). 

For example in the case of 3 Sf 

P ss + = JjJdkoJ d|fe||fc| 2 a; a ^ a (fco,|fc|)| 2 

irr M 1 

= + E *) {(M/2 _ Ek + ier _ kir I*l> 2 • 

If we assume that T 3s +(ko, \ k\) has no pole in the fco-plane (when there is a pole, we 
need careful treatment of the location of the pole of T.), we have the only one pole from 
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the propagator S+, namely at k = 4f — E k + ie. As a result we can obtain 



1 r o f^B S +(k ,k \)\ 

P 3q+ = — / d\k\\k\ 2 — ^ , (A.6) 

35 i N' J 1 11 1 (M 2-E k ) ' 



where ko = M/2 — E k . Here we have used 



2 



i \ „ ,2 t „ a ( i 



a / . i 



27T?-=r 3s +(fco, |fc|) 



d r 2 r w*? 



r 35+ (fco,|fc|) 2 . (A.7) 

4/c 1 



From Eq. (IA.6jl we can consider that 

r 35 +(& , |fc| 



(M/2 - E k 



(Ai 



corresponds to the momentum space wave function in the non-relativistic framework. 
From the Schrodinger equation 

(H + V)ip = Eip , (A.9) 
the wave function can be expressed as 

i> = W ^V^ = -±—T. (A.10) 

Hi — T1Q Hi — I1Q 



Here we have used Vip = T, which is obtained from Eq. f|2.3[) and Eq. (|2.4|) . Comparing 
Eq. (|A.8P and Eq. (jA.lOp . we can immediately understand that Eq. (|A.8P corresponds to 
the Schrodinger wave function. 

Appendix B Negative energy contributions in a 
static potential 

In this appendix, we consider the BS equation in a static separable potential, where 
we study the negative energry contribution for bound and scattering states. 

Let us start with the discussion of the scattering states. To be definite, we consider 
the scattering of the 3 S , J I " channel, where intermediate states are restricted to 3 S± , 3 Df, 
1 Pf and 1 P° , as used in this work. Hence in the BS equation (|2.5p . the in and out states 
are a = (3 = 3 S± . The other channels enters intermediate states labeled by 7 and 5. The 
propagators for the negative energy components are S e and S D which are the combination 

of 5_| and S |_ as given by the Eqs. (|2.7p . The location of the poles of S e and S are 

given in the Fig. [3] in the &o-P iane - 

In the static approximation for the g-function, we set ko = 0: 

lp Z(u in\ W fc =Q W , , 

94 i fc 0, 1*1) - (fc2 _ k2 _ 2 - . e)2 — (k2 + ^ 2)2 , IB.lj 
, 4 (^lkl)= 7 3- (fco2 _ k2 _ /3f+ . e)2 — 0. 
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Substituting these equations, the k integral for the negative energy contributions becomes 

/ dk oj^^{S + Ako,\k\) + S^(k ,\k\)} = 0. (B.2) 

Here, we used the facts that ^-functions are indipendent of ko and the poles of S-\ and 

S (. are located in the same side with respect to the real axis on the complex plane. In 

this way, we can show that the contribution of intermediate states of negative energy do 
not contribute to the scattering amplitude. 

Next, let us consider the bound states. For this purpose we consider the BS equation 
for the bound state amplitude of each channel (^Sf, S D^, Pf, 1 P°)- In particular, we 
need to look at the one of negative energy, say, a = 1 Pf and investigate whether it survives 
or not. Then, the BS equation we should solve is 

0ip«(Po,|p|) = / dk k 2 d\k\ Y Sipf/3(PoJp|) Vg 7 (j> , |p|, feo, |fc|) <Ay(&o, \k\) 

Pi 

N 

= J2 c ihj^kiSi P ep{po, |p|)5 t - (po, IpI) 

/?7<5 i,j,k,l=l 

x I dk k 2 d\k\S lS (k ,\k\)gy\k ,\k\)g[ 5) (k ,\k\). (B.3) 

Here, unlike the BS equation of the form ip = 4> + SVip for scattering states, there is 
no inhomogeneous term correponding to the incoming plane wave <f>. In the second line, 
we use the separable ansatz soulution of Eq. (|2,12p . Using the static approximation for 
the ^-function, the ko integral for the negative energy components becomes zero as in the 
case of the scattering states. However, the part related to the positive energy components 
remains. Therfore, the BS amplitude <f>ipe(po,\p\) does not become zero in the static 
approximation. For a = l P° case, we can discuss in the same way. Therefore, we should 
expand the BS amplitude by the partial wave including the negative energy components. 

From these discussions for scattering and bound states in the static approximation, 
we can justify the prescription of the pole locations as discussed in section (12. 4h . 
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